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We propose a multifractal model for short-term interest rates. The model is a version of the 
Markov-Switching Multifractal (MSM) , which incorporates the well-known level effect observed 
in interest rates. Unlike previously suggested models, the lcvcl-MSM model captures the power- 
law scaling of the structure functions and the slowly decaying dependency in the absolute value 
of returns. We apply the model to the Norwegian Interbank Offered Rate with three months 
maturity (NIBORM3) and the U.S. Treasury Bill with three months maturity (TBM3). The 
performance of the model is compared to level-GARCH models, level-EGARCH models and 
jump-diffusions. For the TBM3 data the multifractal out-performs all the alternatives considered. 



1. Introduction 



Interest rates play an important role for financial institutions, for instance in risk management. Interest- 
rate risk is often assessed via simple stress tests, where one considers parallel shifts in the yield curve 
(typically 100 or 200 basis points) and reports the changes in the value of the portfolio. This approach can 



be improved by stochastic modeling of future movements in the interest-rate yield curve. Vasicek 1977 and 
|Cox et al.| [1985] argue that the yield curve is given by the spot rate alone. Short-term interest rates, such 



as NIBORM3 and TBM3, are frequently used as proxies for the spot rate, and hence accurate modeling of 
these time series is potentially very important. To our knowledge the present paper is the first published 



study considering the Norwegian rate, while the TBM3 has been studied by e.g. Andersen and Lund| [1997 
|Chapman and Pearson| [2001] , [Durham| [2003| , | Johannes| [2004] and |Bali and Wu| [2006| 7~ 



Traditionally, short-term rates have been modeled by Ito stochastic differential equations on the form 



dR( s ) = f(R(s))ds + cR(sydB(s) 



(1) 



where s is the time variable, R(s) denotes the risk-free rate and B(s) is a Brownian motion. If we discretize 
this equation, by letting s = tAs and rt = R(s), then we obtain a stochastic difference equation on the 
form r t = r t -\ + f(r t ~i)At + cAs^-^rJ^wt, where wt are independent and Gaussian distributed random 
variables with unit variance. It is convenient to write this equation on the form 

r t = Ht + &t, (2) 

with fi t = T t -\ + Asf(r t -i) and a t = rj^xt- Here xt = crwt with a = cAs 1//2 . Throughout the paper we 
will consider linear drift terms on the form f(r) = Aq + A\r. If A\ < 0, then ([I]) has a stationary solution, 
and the discretized equation has a stationary solution for sufficiently small As > 0. For a fixed As > we 
write [it = r t -\ + ao + a\r t -\. In this case we have stationarity for — 2 < a.\ < 0. 

The number 7 > is called the Constant Elasticity Variance (CEV) parameter. For 7^0 the models 
feature the so-called level effect. It is generally accepted that this effect is present in interest-rate data, see 



for instance Longstaff et ah, 1992 . The level effect introduces volatility persistence, i.e. strong dependence 
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between the absolute values of the increments Art 
then the variables 

Ar t ~ ("o + am-i) 



r% — rt-i- However, if B{s) is a Brownian motion, 

(3) 



t-i 



are i.i.d. and Gaussian, and hence the volatility persistence will vanish under a simple transformation of 
the increments. Such processes are called pure- level models and have been studied by Cox et al. 19851 and 



Longstaff et al. [1992] . Using time-series data for short-term interest rates we can optimize the likelihood 

6l We can then transform the 
jjc) and[2[c)) are realizations 



to determine the parameters ao, ct\ and 7. These results are shown in table 
data according to (3). The resulting time series (which are plotted in figures 
of the so-called normalized increment process. Just from inspection of figuresnlc) andWc)) we can observe 
that the resulting time series exhibit strong volatility persistence. This is confirmed in figures[TJa) and[T|b) 
where we have plotted the autocorrelation functions for the absolute values of the normalized increments. 
As a consequence of these observations we conclude that pure- level models are insufficient, and we will 
therefore replace the process wt with dependent variables xt- 

In our empirical investigations we find that the dependency structure in the variables xt resembles the 
stylized facts of logarithmic returns of asset prices, namely that the variables themselves are uncorrelated 
(or weakly correlated) whereas their absolute values have slowly decaying autocorrelation functions. To 
describe this dependency several authors (e.g. Brenner et al. 1996] and Koedijk et al. [1997] ) have suggested 
using a generalized autoregressive conditional heteroskedasticity (GARCH) model for the process xt- The 
corresponding process rt is then called a level-GARCH model. These models can be further improved by 
letting x t be an Exponential GARCH (EGARCH) model or a Markov Switchin g GARCH (MSGARCH) 
model. The level-EGARCH model was introduced by Andersen and Lund [1997 and has been reported to 



perform better than the standard GARCH model on the TBM3 data. 

In this work we propose a multifractal model, specifically a level-MSM model, as an alternative to 
level-EGARCH models and level-MSGARCH models for short-term interest rates. This model is a slight 
modification of the standard MSM model constructed by Calvet and Fisher (2004). The motivation for 
introducing multifractal models for short-term interest rates is similar to the motivation for multifractal 
modeling of asset prices, namely that these models capture the dependency structure and scaling properties 
of the process xt- Secondly, encouraging results (in- and out-of-sample) have been obtained in a preliminary 
study of the NIBORM3 data [L0vsletten[ |Autumn 2010] 



We follow Calvet and Fisher (2004) and use an adjusted version of the Vuong test Vuong, 1989 for 
model selection, and the main result of this work is that the level-MSM model performs well compared 
against level-EGARCH models, level-MSGARCH models and jump-diffusions. The paper is organized as 
follows: In section [2] we define the level-MSM model and consider some basic properties. Technical details 
about multifractal processes are presented in appendix [I] A brief description of the alternative models are 
presented in section [3j and in section [4] the results of the Vuong test are presented. In section [5] we make 
some concluding remarks. 



Remark 1.1. The data analyzed in this paper is freely available online. Both data sets are given with daily 
resolution. The TBM3 data is taken from the period 1951-01-04 to 2010-09-22, and consists of 14172 data 
points. The NIBORM3 data is taken from the period 1986-01-02 to 2010-09-24, and consists of 6231 data 
points. The TBM3 data contains several days with zero value, which causes certain technical problems in the 
stochastic modeling. The problem is resolved by a simple variable shift r 1— > r + b. In order to be consistent 
with our set of models, we estimate the constant b > by estimating the conditional standard deviations 
sdv[Art|rt_i = r] (for various r ) and fitting a function cr 1 — b to this data set. With this approach we find 
6 = 0.03. 



2. Multifractal models 



The application of multifractal processes to finance was introduced by Mandelbrot et al. 1997 . A process 
X(t) is multifractal if the structure functions S q (t) = ¥,[\X(t)\ q ] are power laws in t, and one define a 
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scaling function £(g) through the relation S q (t) ~ as t — )■ 0, i.e. 

sw t^O log* 

Using Holder's inequality it is easy to show that the scaling function £(g) is concave, and we also note that 

if the process X(t) is h-sett similar, i.e. X(at) = a h X(t), then ((q) = hq. We are therefore interested in the 
situations where ((q) is strictly concave. In this case we see that 

E[X(t)*] 



E[X(t) 



212 



t «4)-2C(2) ^ ag t _ >0 _ 



If we assume that X(t) also has stationary increments, then this implies that the At-lagged increments 
X(t + At) — X(t) are more leptokurtic for small At than for larger At. In particular, X(t) can not be a 
Gaussian process. 

In financial time series one will most often want the process X(t) to have uncorrelated and stationary 
increments. In order for this to be satisfied we must impose the condition ("(2) = 1, because otherwise the 
variables AX{t) = X(t + 1) — X(l) have slowly decaying autocorrelation. In fact, for £(2) ^ 1 

Cov(AX(0), AX(t)) ~ t 2Hd - 2 , 

where := £(2)/2. The advantage of a multifractal process (with strictly concave scaling function) is 
that it has strongly dependent increments, even for Hd = 1/2. For instance, in the MSM model 1 we have 
C(2) = 1 and 

Cov(|AX(0)| 9 , \AX{t)\ q ) ~ t-KWi-WM)) . (4) 



See proposition 1 in |Calvet and Fisher, 20041 . This inherent long-range volatility persistence serves as our 



motivation for modeling short-term interest rates using multifractals. 
2.1. Mandelbrot's MMAR processes 

The multifractal models introduced by Mandelbrot et ah] [1997| are stochastic processes on the form 



X(t) = aT 1 ' 2 B{Q(t)), (5) 

defined for < t < T, where G(t) = /i([0,t]) is a random probability measure on [0,T] and B(t) is a 
Brownian motion with E[B(1) 2 ] = 1. The process G(t) is itself a multifractal process with scaling function 
(o(q), and if \x and B(-) are independent, then Cx(l) = By construction, the measures we will 

mention in this paper satisfy E[/^([0,t])] oc t, i.e. Ce(l) = 1- It follows that Of (2) = 1. We note that 
processes with correlated increments are useful in other applications (see e.g. |Rypdal and Rypdal 2010| ), 



and the Borownian motion B(t) can in these situations be replaces by a fractional Brownian motion with 
H^l/2. 

There are several choices for the measure fi, and for completeness we give three examples in the 
appendix. We also show how a particular construction, the Poisson multifractal, leads to a discrete-time 
MSM model defined in the next section. 

2.2. The MSM model 

The Markov Switching Multifractal model of order K > 1 is given by 

x t = ag(M t ) 1 / 2 e t , (6) 



1 It is pointed out by Lux 2006 that the MSM model is only a finite- level approximation, an d so equation (jjtjl is only valid on 
time scales up to b K . Here b > and K 6 N are parameters in the MSM model. See section 2.2 
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where Mj = [M\$, ■ ■ ■ , Mftr,*)) and the function g : R K ^ R is the product of the vector components. 
The innovatons £t are Gaussian distributed with unit variance, and in this setting also assumed to be 
independent. The process Mt is a Markov chain defined by the following updating scheme: 

J drawn from M with probability 
Mk t = s (7) 
I Mk t t-i with probability 1 — A&. 

Components are updated independently of all previous updates, and the frequencies A& are related to each 



other through (A.4). 



For simplicity we choose a two-point distribution M: 

P(M = m ) = P(M = mi) = - , with m G (1, 2) and m x = 2 - m . (8) 



This version of the MSM model is known as the binomial multifractal |Calvet and Fisher, 2004 . To higlight 



the number of multipliers K we will use the notation MSM(X). Note that, unlike GARCH(j>, q), the number 
of parameters remains unchanged with increasing order. 

Let {m ! | i = 1, . . . , d := 2 K } denote the sample space of M^. We also define vectors lo(x) : M — > M. d 
and pt : N — > M. d with components C0k( x ) = n(x\a 2 g(m k )) and pkj = 1P(Mj = m fc |xi) respectivly. Here 
x t = (x2, ■ ■ ■ ,xt), and n(-\a) is the density of the normal distribution with zero mean and variance a 2 . 
With this notation the density of x t |x t _i is given by 

p(x t \x.t-\) = pJ-i Au; (xt), (9) 

where A is the transitition matrix of the Markov chain, that is Aj k = P(Mt = m fe |Mi_i = mP). This is 
seen by conditioning on the underlying Markov-chain at time t — 1 and t, and using Bayes' rule. It is easily 
seen, once again using Bayes' rule, that the vectors pt follow the recursion 

pj = — ^ T (^) * (p* t -i^), (io) 

p(x t \xt-i) 

where * denotes the Hadamard product. We start the recursion with the limiting probabilities of Mj. 

2.3. The level-MSM model for interest rates 

We propose the following level-MSM model for short-term interest rates: 

n -n-i = ao + air t -i +r]_ 1 x t , (11) 
where xt is the binomial MSM model defined above. The likelihood 

n 

L = 2^1og/(r t |r t _i) , 

t=2 

for data ?"2, . . . ,r n , with r\ taken as a pre-sample value, now follows from equations ^ and (10) together 
with the relation 

/(r t |r t _i) = -j— p(x t \x t -i) , 
r t-i 

where / is the density of rt|r t _i. 

We have fitted the level-MSM models to the two time series under consideration for K = 2, . . . , 9. For 
the case K = 9, the maximum likelihood (ML) estimates for the parameters ao and a\ are reported in table 
|6| In this table the estimates of ao and a\ for the alternative models are also included. As expected the 
estimates for the parameter a\ are negative, but we observe that none of these ML estimates are significantly 
different from zero for the TBM3 data. The estimates for ol\ are also very small for the NIBORM3 data, 



and we will therefore consider the model defined by (11) with ot\ = 0. The corresponding ML estimates for 



the other parameters are presented in table [6j From these results we also see that the likelihoods increase 
monotonically with K. Using a Vuong test we compare the level-MSM of order K = 9 against level- 
MSM models of lower order. When comparing two models with this method, the null-hypothesis is that 
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both models are equally far from the data-generating process measured by the Kullback-Leibler distance. 
Hence small p- values indicate that the level-MSM model of order K = 9 is significantly closer to the true 
generating process than the level-MSM models of order K = 7 or lower. 

The same model-selection-test is used to compare the level-MSM model with various alternative models, 
and more details on the test are presented in section [4} 



3. Alternative models 

In this section we briefly discuss some processes which we use for benchmarking the level-MSM model. 
As for the level-MSM model we have first considered these models with drift terms on the form fit = 
rt-i + ao + a\rt-\. The ML estimates for the parameters ao and a\ are reported in table [6} Again we see 
that the contribution from the parameter a\ is negligible. Consequently we will consider models with drift 
terms on the form fit = Tt—\ + ao- 



3.1. The level- G ARCH model 



In this paper we prefer the standard level-GARCH model proposed by Koedijk et al. 11997 . As is common 



for this model, we include student-i distributed innovations. The model then reads 

n = ao + n_i + cr t , G t = h^ 2 fl_ x et, h = oq + a\xl_ x + bht-i, (12) 

where the innovations et are i.i. d. i„(0, l) 2 . For 7 = we have the pure-GARCH(l,l) model, and for 
a\ = b = we have the standard level model. 

The estimated parameters for the NIBORM3 and TBM3 data are presented in table [6j We observe that 
simulated paths with these exponents have far too wild fluctuations compared to the real data, indicating 
that this model fails to accurately describe the interest-rate fluctuations. In addition, we know that the 
GARCH models exhibit exponentially decaying autocorrelation functions for the absolute values of the 
increments. This means that the long-range volatility persistence observed in the short-term interest rates 
is not inherent in these models. 



3.2. The level-EGARCH model 

Instead of using (12), Andersen and Lund [1997 propose using the EGARCH model. They find that this 
model gives an adequate fit to the TBM3 data, and a better fit compared to the standard GARCH model. 
In the EGARCH model the logarithm of the conditional variance replaces the conditional variance. The 
variance-recursion is then 

log h t = a + ai£ t -i + o 2 |E t _i| + &log/t t -i- (13) 

The extra parameter a\ controls potentially different responses to positive and negative returns. The use of 
the logarithm guarantees positive values of the volatility for all parameter values. In addition, this model 
provides the extra flexibility by letting the conditional distribution be non-symmetric. As we will become 
apparent from the results presented in the next section, the EGARCH model gives better results than the 
GARCH model for both of the time series considered in this paper. This confirms the results of Andersen 
and Lund|p97| . 

The ML estimates for the parameters in the EGARCH model are presented in table [6| 



3.3. Jump-diffusions 

The final class of alternative models considered are the jump- diffusions. Both Johannes 



2004 and Das 



[2002] propose these processes in order to describe the large spikes observed in interest-rate data. As a 
discretized jump-diffusion model we use the following specification: 



n 



a + r t -i + r7_i(a?t + J t z t ) 



2 i;,(0, 1) is a centralized, unit- variance student-f distribution with v degrees of freedom. 



() 



where Jt is the jump-indicator assumed to follow a Bernoulli-distribution. The probability of a jump taking 



place at time t is then given by (1 + exp(— c — dr t -i)) . By letting xt = h 



1/2 



conditional heteroschedacity in levels is also accomodated Das , 2002 Hong et al. , 2004 



£t follow a GARCH-process, 
The innovations 



are distributed as zt ~ AA(0,r 2 ) and et ~ A/"(0, l).The variance recursion in the GARCH process is given 
by 

/ Ar t _i -a V . , , 
ht = ao + ai ~ +bh t -i, 



\ 't-2 / 

where Art = r t — r t-i- 

The ML estimates for the parameters in the jump-diffusion model are presented in table [6j 
4. In-sample comparision 

To test the binomial MSM model against the alternative models we emply a version of the Vuong test 



which is adjusted for heteroschedacity and autocorrelation (HAC) |Calvet and Fisher 2004| . For each of 
the alternative models the null-hypothesis is that this model and the level-MSM model of order K = 9 are 
equidistant from the true data-generating process, measured in the Kullback-Leibler distance |Kullback 



and Leibler, 1951 . In the classical Vuong test it is assumed that the data-generating process is i.i.d., and 
then the log-ratio of the likelihoods for the two models will converge to a normal distribution with zero 
mean, and with a variance which is consistently estimated by the sample variance for the log-ratio of the 
likelihoods. Using the corresponding normal distribution one can then easily calculate a p- value under 
the null- hypothesis. In the HAC-adjusted version of this test, the data-generating process may exhibit 
dependence, but the variables rt should be identically distributed. This is satisfied if we assume that the 
data-generating process rt is stationary. 

The results of this test are presented in table |6| We observe that for the TBM3 data, the level-MSM 
model of order K = 9 is significantly closer to the data-generating process than any of the alternative mod- 
els. For the NIBORM3 data the EGARCH model performs best, whereas the level-MSM model performs 
better than the standard level-GARCH model and jump-diffusions. 



5. Concluding remarks 

In this paper we have introduced a multifractal model for short-term interest rates. The model combines the 



well-established level effect described in Longstaff et al. , 1992] with the discretized multifractal model of 



Calvet and Fisher, 2001 . In a comparison with level-GARCH, level-EGARCH and jump-diffusions, we find 
that this model well describes the fluctuations of the TBM3 and NIBORM3 time series. The main result 
of this work is that the level-MSM outperforms all alternatives for the TBM3 data. This motivates further 
research on multifractal modeling of short-term interest rates, in particular an out-of-sample analysis of 
the level-MSM model. 

It is also interesting to note in the level-GARCH model the parameter estimates (for both the TBM3 
data and the NIBORM3 data) fall outside the covariance-stationarity region. As a result the level-GARCH 
model has a wild volatility pattern, which does not seem to be an accurate description of the interest-rate 



data. This confirms the results of e.g. Andersen and Lund 1997 



Appendix A Examples of multifractal measures 



In this appendix we give examples of random measures that can be used to define multifractal processes 
via ([5]). For the Poisson multifractal we show how a discretization leads to the model in ([6]). 

Example 1: In the simplest case fi is a randomized dyadic Bernoulli measure with probabilities pi = p 
and P2 = 1 — P- This measure is constructed through an iterative procedure, where we in the first step 
divide the interval [0, T] in two pieces Aq and Ai of equal length. One of the intervals is chosen at random 
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and given a mass p% fj,([0,T]), while the other interval is given the mass p 2 /i([0,T]). This procedure is then 
repeated recursively, i.e. the interval Ao is divided into the equally sized intervals Aoo and Aoi- One of 
these intervals is given the mass pf fi([0,T]) and the other is given the mass P1P2 /-t([0, T]). Formally the 
measure can be defined by letting / be a random bijection 3 {0, 1} — > {pi,P2} and fi lt ... : i k be independent 
copies of/. The measure is constructed by assigning the mass fi(^i u i^...,i k ) = /{iijfhfa) ■ ■ • fi 1 ,i3,...,i h _ 1 (ik) 
to the dyadic intervals 

Ai!,...,^ = [T • (0.*i • • • i k ) 2 , T ■ (Oil • • • i*)a + T • 2~ fc ] , i n G {0, 1} . 

A simple combinatory argument shows that for t = T ■ 2~ k , the random variable 0(t) has density 

k 



= i E (f) *(* " P n ( x " • 



n=0 

When the processes B(t) and ®(t) are independent, we obtain the density of X(t): 



Px,t( x ) 



/ P CT Ti/2 B , s (^)P0,t(s)^ = iE 



^ fk\ I 2<7 2 (Tp' l (l-p) fc -") 



n=o vv J2vra 2 (Tp"(l - p) k ~ n ) 



From this density we can easily calculate the structure functions S q (t) = ~E[\X(t)\ 9 ] and see that S q (t) ~ t^ q > 
ast^O, where ((q) = 1 - T(g/2) and T(q) = \og 2 (p q + (1 - p) 9 ). 

Remark 5.1. The function T(q) is sometimes called the scaling function of the random measure. The dyadic 
measure defined by /i(A^ 1 . i fc ) = Pi ± - ■ ■ Pi k has a multifractal spectrum f{a) given by the Legendre transform 
ofT(q): f(a) = mi q ^ {aq — T(q)}, and the Hentchel-Procaccia dimension spectrum is D q = (q — l)~ 1 T(q). 
See e.g. fPesin , 199^ for an account of the relation between scaling functions and spectra of fractal dimen- 



sions. The function a h-> /(a/2) is often called the singularity spectrum of the process X(t). 
Example 2: A different class of multifractal measures are the 6-adic random multiplicative cascades (see 



e.g. Mandelbrot et al. , 1997 for a more detailed account). For an integer b > 2 and i n E {0, 1, . . . , b — 1} 
we define the 6-adic subintervals of [0, T] as 

A n ,..., ifc = [T ■ (Oil . . . i k ) b , T • (O.ii . . . i k ) b + T ■ b- k ] . 

For a positive random variable M, with E[M] = let 

M\Ajj ) — Mj-^Mi-^^ ' ' ' -^ii«2"-ifc ^il,.-.,ifc ) 

where Mj 1) ... ) j n are independent copies of M and 

^u,...,i fe = lh 11 'y ] Mi x • • • Mi 1 i 2 ...i k Mi 1 i 2 ...i k j 1 • • • Mi 1 i 2 ...i k j k ...j n . 

jl,-,jn 

Again a stochastic process is constructed according to equation ^ with @(t) = £i([0, t]). This process 
has a scaling function ("(g) = 1 — T(q), where T(q) = log&E[M 9 ]. Popular choices for the multiplier M are 
log- normal distributions, which give quadratic scaling functions, or other log-infinitely divisible distribu- 
tions. Almost every realization of the measure /U has the multifracal spectrum f(a) = mf q ^{aq — T(q)}. 

Example 3: The Poisson multifractal measure on [0, T] generalizes the 6-adic multiplicative cascade by 
introducing randomness in the construction of the intervals A^ j 2 ... In the original multifractal models, 
each interval on level k are divided into b pieces of equal length at level k + 1. As a result the interval 
length decreases as b~ k with the level k. In the Poisson multifractal measure, the splitting of an interval at 
level k is preformed by drawing the lengths of the new pieces randomly from an exponential distribution 
with rate l k +i- This means that the mean and median length of an interval at level k + 1 are 1/Zfc+i and 



3 With probability 1/2 for each of the two outcomes (0, 1) i-> (pi,P2) and (0, 1) i-> (f»2,Pl) 
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log(2)/Zfc + i respectively, so to maintain exponential decay of interval lengths (as a function of level), one 
chooses If. = b k ~ x l\. Note that b no longer is restricted to the integers. 

Formally, the Poisson multifractal measure on the interval [0, T] is defined via a sequence of measures 
Hk specified on randomly generated intervals 



Aii,i2,...,t* — [iii,i2,... ,iki ^ti,i2i—>(*fc+ 1 )] ' 
i k are defined by the following recursive construction: Let {r^ 



where the numbers 
independent and exponentially distributed random variables with rates l k 



b k Define 



Ni- 



i*fe-i 



max < m : T h , 



3=1 



•i*fc- 



_ x (j) < diam(Aj 1) » 



-■)} 



denote 



(A.l) 



and 



til , 



•i»fc 



>*fc-i 



j = l T *l, 



if 
if 



(A.2) 



ii,...,(i fc _i+l) 



This means that the interval A,- 



. i k is divided into 



process with rate l k +i- We start with an interval A 
via the formula 



in = 

1 < i fc < AT il> ... )ifc _ 1 . 
if i fc = JV^,...,^ + 1 

subintervals by the cuts made by a Poisson 
0, T]. A sequence fi k of measures can now be defined 



T 1 diam(Aj li 



) M h M iui2 ■ ■ ■ M iu ..„ 



1-h 1 



where M^...^ are independent copies of a positive random variable M satisfying E[M] = 1. If E[M 2 ] < 6, 
then /ifc converges weakly to a Borel measure /i on [0, T\. With this choice of random measure the model 
given by equation ^ has scaling function 

C(g) = 1 - T(q/2) , where T(q) = 1 - q + log fe E[M q ] . (A.3) 

Calvet and Fisher constructed their MSM model by discretizing the time interval [0, T] and assigning 
discrete geometric (rather than an exponential) distributions on the waiting times r. More precisely one 
will fix an integer K > 1 determining the number of levels that are to be included in the discrete model, 



and consider the integer values {0,1,..., T}, where T 



m 



K 



for some positive integer m. As for the 



construction of the Poisson multifractal one makes random partitions Aj x 



[tii, 



of the interval [0, T\. These are constructed using the recursive procedure described by equations (A.l) 

i k are replaced by discrete random variables fj li ... i j fc with geometric 

T / T . This choice of parameters 



and (A.2), where the variables 

r') 



distributions F(r iu ... )ih = r') = X k (1 - A fc ) r ' _1 , with 1 - X k = e^ 6 ' 
implies that the median of fi lt ,,,^ k is proportional to ^ 1 6 _n . Note that the parameters are given by Ai 
through the formula 

X k = 1 - (1 - AO 6 "" 1 . (A.4) 

A measure jj, is defined by specifying that if t is an integer and [t — l,t] C Ai u ...,i K , then 

- 1, t]) = m- K M h M lui2 ■ ■ ■ M ilt ^ K , 

where Mi ± i k are independent copies of a positive random variable M with E[M] = 1. Again we define 
Q(t) = p,([0,t]), and by taking the composition with a fractional Brownian motion we get a discrete-time 
stochastic process X(t) = C B(@(t)). We denote xt = X(t) — X(t — 1), and observe that 

x t ^Cil{[t,t-l]) l l 2 (B(t)-B(t 



1 



a (M h M iui2 ■ ■ ■ M ht 



where a = C m K and e< = B{t) — B(t — 1) is a discrete version of a white Gaussian noise. We can simplify 
notations by denoting M ky t = Mi u ... t i k for t G Aj^...^. 

We finally remark that there exists a large class of multifractal random measures, knowns as log- 
infinitely divisible cascades Bacry et al. 2001 Bacry and Muzy, 2003; Bacry et al. 2008 , which have 
multifractal scaling and stationary increments. However, there only exists approximate maximum likelihood 



<) 



methods for these processes [L0vsletten and Rypdal , 2011] , and hence the Voung-testing preformed in this 
paper is not available. 



Acknowledgment. This project was partly funded by Sparebank 1 Nord-Norge and the Norwegian Re- 
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6. Figures and tables 
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Fig. 1. (a): The NIBORM3 data r t for the time period 1986-01-02 to 2010-09-24. (b): The one-day increments Ar t = r t -r t -i- 
(c): The normalized one-day increments Art/r^i- The value of the CEV parameter is 7 = 1.61. 
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Fig. 2. (a): The TBM3 data r t for the time period 1954-01-04 to 2010-09-22. (b): The one-day increments Art = n - r t -i- 
(c): The normalized one-day increments Art/rJ_ 1 . The value of the CEV parameter is 7 = 0.38. 
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for NIBORM3 
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Fig. 3. (a) : Autocorrelation functions for one-day increments in the NIB0RM3 data. The variables considered are the standard 
increments Art = ft ~ ft— 1, the normalized increments Art/r7i an d the absolute values of these two time series. The value 
of the CEV parameter is 7 = 1.61. (b): Same as (a), but now for the TBM3 data. The value of the CEV parameter is 7 = 0.38. 
In both figure the dotted lines represent a 0.95 confidence interval for the autocorrelation function assuming independence. 
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Innovations 



Table 1. CEV-models. 
Qg 7 a 



2.01 



logi 



U.S. Treasury Bill 
Normal 0.0009709 0.3755 0.04788 

t u -0.0002434 0.5377 0.2186 

NIBORM3 

Normal 0.0002863 1.61 0.006298 

t u 0.0002419 0.988 0.0921 



15707.94 
2.01 20688.54 



4771.78 
8041.42 



Note: The constant elasticity volatility (CEV) model of [Longstaff 
et al. 1992 with normal and student-t innovations. 



K 10 J x a 



;»0 



Table 2. Multifractal models 
b Ak a logL 



BIC 



Vuong HAC-adj. 



U.S. Treasury Bill M3 



2 


0.2752 


0.2824 


1.808 


23.85 


0.1391 


0.07037 


22309.49 


-3.1446 


13.573 


9.091 




(0.247) 


(0.0106) 


(0.00471) 


(4.92) 


(0.011) 


(0.00162) 






(< 0.001) 


(< 0.001) 


3 


0.2052 


0.1913 


1.73 


10.62 


0.1660 


0.08078 


22676.53 


-3.1964 


11.618 


10.138 




(0.234) 


(0.0175) 


(0.0062) 


(1.45) 


(0.0163) 


(0.00336) 






(< 0.001) 


(< 0.001) 


4 


0.2276 


0.2269 


1.661 


9.643 


0.2956 


0.07202 


22818.9 


-3.2165 


8.4076 


7.562 




(0.226) 


(0.032) 


(0.00674) 


(1.23) 


(0.0462) 


(0.00357) 






(< 0.001) 


(< 0.001) 


5 


0.0643 


0.2280 


1.609 


9.804 


0.6761 


0.06865 


22901.03 


-3.2281 


5.036 


4.805 




(0.222) 


(0.0191) 


(0.00667) 


(1.04) 


(0.0597) 


(0.00265) 






(< 0.001) 


(< 0.001) 


6 


0.0375 


0.2115 


1.559 


6.577 


0.778 


0.06368 


22920.57 


-3.2308 


4.346 


4.335 




(0.219) 


(0.0282) 


(0.00744) 


(0.564) 


(0.0541) 


(0.00360) 






(< 0.001) 


(< 0.001) 


7 


0.0660 


0.1744 


1.527 


5.713 


0.7928 


0.08208 


22933.55 


-3.2326 


3.4883 


3.112 




(0.223) 


(0.0291) 


(0.00752) 


(0.465) 


(0.0561) 


(0.00441) 






(< 0.001) 


(0.001) 


8 


0.0839 


0.1575 


1.503 


4.932 


0.8474 


0.06956 


22951.89 


-3.2352 


1.131 


0.938 




(0.215) 


(0.0260) 


(0.00815) 


(0.355) 


(0.0509) 


(0.0036) 






(0.129) 


(0.174) 


9 


0.0693 


0.1984 


1.462 


3.864 


0.931 


0.06029 


22955.80 


-3.2358 








(0.215) 


(0.037) 


(0.00724) 


(0.255) 


(0.0396) 


(0.00443) 




















NIBORM3 










2 


1.156 


1.805 


1.843 


16.23 


0.1016 


0.006092 


8514.45 


-2.7250 


5.219 


5.298 




(0.408) 


(0.04) 


(0.00547) 


(4.68) 


(0.0116) 


(0.000516) 






(< 0.001) 


(< 0.001) 


3 


1.041 


1.761 


1.772 


18.37 


0.4259 


0.005674 


8671.87 


-2.7755 


3.054 


3.093 




(0.378) 


(0.0391) 


(0.00685) 


(3.07) 


(0.0459) 


(0.000463) 






(0.001) 


(< 0.001) 


4 


1.020 


1.482 


1.717 


12.43 


0.3939 


0.01401 


8760.38 


-2.8039 


3.424 


3.588 




(0.378) 


(0.043) 


(0.00794) 


(1.93) 


(0.0483) 


(0.00129) 






(< 0.001) 


(< 0.001) 


5 


1.13 


1.933 


1.642 


7.495 


0.5455 


0.004534 


8779.88 


-2.8102 


3.292 


3.265 




(0.381) 


(0.0635) 


(0.0100) 


(1.04) 


(0.086) 


(0.000555) 






(< 0.001) 


(0.001) 


6 


1.079 


1.620 


1.590 


5.347 


0.6584 


0.00897 


8789.49 


-2.8133 


2.702 


2.866 




(0.379) 


(0.0481) 


(0.0107) 


(0.542) 


(0.0726) 


(0.000841) 






(0.003) 


(0.002) 


7 


1.006 


1.934 


1.537 


4.431 


0.9000 


0.003719 


8798.81 


-2.8162 


1.809 


1.727 




(0.379) 


(0.0685) 


(0.00926) 


(0.448) 


(0.0748) 


(0.000464) 






(0.0352) 


(0.042) 


8 


1.023 


1.985 


1.500 


3.536 


0.9173 


0.003220 


8802.66 


-2.8175 


1.058 


0.940 




(0.380) 


(0.063) 


(0.00933) 


(0.292) 


(0.0624) 


(0.000388) 






(0.145) 


(0.173) 


9 


0.998 


1.882 


1.464 


2.964 


0.949 


0.003706 


8804.44 


-2.8181 








(0.379) 


(0.112) 


(0.00927) 


(0.216) 


(0.0484) 


(0.00102) 











Note: ML-estimates for the binomial multifractal with estimated standard errors in parantesis. The Vuong-column reports 



the likelihood ratio statistic with corresponding p- value in brackets Vuong 1989| . The null hypothesis is that MSM(9) 
and MSM(K) have equal Bayesian Information Criteria, with the alternative hypothesis being that the MSM(9) is closer 
to the true data generating process. The HAC-adjusted version of the Vuong test jCalvet and Fisher] |2004] corrects for 
heteroschedacity and autocorrelation in the addends. 



Table 3. Level-GARCH 



10 3 X a 


7 


10 5 x a 


a\ b 


V 


logL 


BIC 






U.S. 


. Treausury Bill M3 








0.0904 


0.1699 


0.7092 


0.1301 0.8915 


3.7995 


22873.03 


-3.2241 


(0.225) 


(0.0496) 


(0.139) 


(0.0092) (0.0062) 


(0.134) 












NIBORM3 








0.9104 


0.9717 


0.4889 


0.20993 0.80838 


3.3305 


8783.27 


-2.8113 


(0.397) 


(0.0719) 


(0.158) 


(0.0209) (0.0135) 


(0.155) 







Note: ML estimates for the parameters in the level-GARCH model. Standard deviations are 
in brackets. 



Table 4. Level-EGARCH 



10 3 X a 


7 


a 


a\ a 2 


b 


V 


logL 


BIC 








U.S. Treausury Bill M3 








0.1039 


0.3441 


-0.2360 


-0.03127 0.2322 


0.9888 


3.866 


22936.53 


-3.2324 


(0.212) 


(0.0648) 


(0.0165) 


(0.00634) (0.0123) 


(0.00154) 


(0.138) 






1.120 
(0.396) 


0.9101 
(0.0887) 


-0.4024 
(0.0372) 


NIBORM3 
0.02930 0.2664 
(0.0104) (0.0188) 


0.9741 
(0.00333) 


3.358 
(0.156) 


8828.90 


-2.8245 



Note: ML estimates for the parameters in the level-EGARCH model. Standard deviations are in brackets. 



Table 5. Jump-diffusion 



10 3 X a 


7 


10 6 x a 


Ol 


b c 


d 


T 


logL 


BIC 










U.S. Treasury Bill M3 










0.464 


0.1670 


5.131 


0.0931 


0.8773 -3.751 


0.2761 


0.08622 


22597.89 


-3.1839 


(0.232) 


(0.0374) 


(0.867) 


(0.0055) 


(0.0062) (0.170) 


(0.0271) 


(0.00624) 






1.085 
(0.417) 


1.872 
(0.0892) 


0.0436 
(0.0198) 


0.1291 
(0.0099) 


NIBORM3 
0.8301 -1.883 
(0.0111) (0.328) 


-0.2119 
(0.0417) 


0.01127 
(0.00216) 


8579.61 


-2.7431 



Note: ML estimates for the parameters in the jump-diffusion model. Standard deviations are in brackets. 



Table 6. In-sample model comparison. 
Model dim(0) logL BIC Vuong Hac-adj. 







U.S. 


Treasury Bill M3 




MSM(9) 


6 


22955.80 


-3.2358 






GARCH 


6 


22873.03 


-3.2241 


-4.064(< 0.001) 


-3.477(< 0.001) 


EGARCH 


7 


22936.53 


-3.2324 


-1.420 (0.078) 


-1.361 (0.087) 


Jump-diffusion 


7 


22597.89 


-3.1839 - 


9.878 (< 0.001) 


-9.240 (< 0.001) 








NIBORM3 






MSM(9) 


6 


8804.44 


-2.8181 






GARCH 


6 


8783.27 


-2.8113 


-1.099(0.136) 


-1.086(0.139) 


EGARCH 


7 


8828.90 


-2.8245 


1.231(0.891 ) 


1.470 (0.929) 


Jump-diffusion 


7 


8579.61 


-2.7430 


-6.964(< 0.001) 


-5.765(< 0.001) 



Note: The Vuong-column reports the test-statistic for differences in BIC. The null hypothesis 
is that the multifractal and the alternative model are equally good, with the alternative 
hypothesis being that the multifractal model performs best. The HAC-adjusted column 
adjusts for heteroschedacity and autocorrelation in the addends. Corresponding p-values 
are in brackets. 



Table 7. Linear drift 



Model 


10 3 x a 


10 4 X Ql 


U.S 


. Treasury Bill M3 




MSM(9) 


0.1214(0.354) - 


-0.1905(1.04) 


GARCH 


0.2821 (0.386) - 


-0.6731(1.10) 


EGARCH 


0.3595(0.310) - 


-1.3524(0.92) 


Jump-diffusion 


0.7074(0.397) - 


-0.8695(1.15) 




NIBORM3 




MSM(9) 


2.228(0.726) 


-3.201 (1.60) 


GARCH 


2.424 (0.746) 


-3.986(1.66) 


EGARCH 


2.515(0.430) 


-3.626 (0.87) 


Jump- diffusion 


2.621 (0.788) 


-4.053(1.77) 



Note: A linear term in the drift term was added and 
the parameters estimated using ML. In the TBM3 
none of the parameters ot\ are significantly differ- 
ent from zero. Standard deviations are reported in 
brackets. 



16 REFERENCES 



References 

T. G. Andersen and J. Lund. Estimating continuous-time stochastic volatility models of the short-term 

interest rate. Journal of Econometrics, 77(2) :343 - 377, 1997. 
E. Bacry and J. F. Muzy. Log-Infinitely Divisible Multifractal Processes. Communications in Mathematical 

Physics, 236(3):449-475, 2003. 
E. Bacry, J. Delour, and J. F. Muzy. Multifractal random walk. Physical Review E, 2001. 

E. Bacry, A Kozhemyak, and J. F. Muzy. Continuous cascade models for asset returns. Journal of Economic 

Dynamics and Control, 32(1):156-199, 2008. 
T. G. Bali and L. Wu. A comprehensive analysis of the short-term interest-rate dynamics. Journal of 

Banking and Finance, 30(4): 1269 - 1290, 2006. 
R. J. Brenner, R. H. Harjes, and K. F. Kroner. Another look at models of the short-term interest rate. 

Journal of Financial and Quantitative Analysis, 31:85-107, 1996. 
L. Calvet and A. Fisher. Forecasting multifractal volatility. Journal of Econometrics, 105(1):27 - 58, 2001. 
L. Calvet and A. Fisher. How to Forecast Long-Run Volatility: Regime Switching and the Estimation of 

Multifractal Processes. Journal of Financial Econometrics, 2(1):49— 83, 2004. 
D. A. Chapman and N. D. Pearson. Recent advances in estimating term-structure models. Financial 

Analysts Journal, 57:77-95, 2001. 
J. C. Cox, J. E. Ingersoll Jr., and S. A. Ross. A theory of the term structure of interest rates. Econometrica, 

53:385-407, 1985. 

S. R. Das. The surprise element: jumps in interest rates. Journal of Econometrics, 106(1):27 - 65, 2002. 
G. B. Durham. Likelihood-based specification analysis of continuous-time models of the short-term interest 

rate. Journal of Financial Economics, 70(3):463 - 487, 2003. 
Y. Hong, H. Li, and F. Zhao. Out-of-sample performance of discrete-time spot interest rate models. Journal 

of Business and Economic Statistics, 22(4):457-473, 2004. 
M. Johannes. The statistical and economic role of jumps in continuous-time interest rate models. The 

Journal of Finance, 59(l):227-260, 2004. 
K. G. Koedijk, J. A. Nissen, P. C. Schotman, and C. P. Wolff. The Dynamics of Short-Term Interest Rate 

Volatility Reconsidered. European Finance Review, 1(1):105-130, 1997. 
S. Kullback and R. A. Leibler. On information and sufficiency. Ann. Math. Stat., 22:79-86, 1951. 

F. A. Longstaff, K. C. Chan, G. Andrew Karolyi, and A. B. Sanders. An empirical comparison of alternative 

models of the short-term interest rate. The Journal of Finance, 47:1209-1227, 1992. 
O. L0vsletten and M. Rypdal. Approximate Maximum Likelihood Estimation in Multifractal Random 
Walks, preprint, 2011. 

O. L0vsletten. Empirical analysis and stochastic modeling of temporal fluctuations in the norwe- 

gian interbank offered rate. Master thesis, written in Norwegian, Autumn 2010. URL [httpl 

/ /complexityandplasmas.net/nordforsk/Papers.html, 
T. Lux. The markov-switching multifractal model of asset returns : Gmm estimation and linear forecasting 

of volatility. Economics Working Papers, 17, 2006. 
B. Mandelbrot, A. Fisher, and L. Calvet. A Multifractal Model of Asset Returns. Cowles Foundation for 

Research in Economics, 1997. 
Y. Pesin. Dimension Theory in Dynamical Systems. Chicago Lectures in Mathematics, 1997. 
M. Rypdal and K. Rypdal. Stochastic modeling of the ae index and its relation to fluctuations in bz of the 

imf on time scales shorter than substorm duration. Journal of Geophysical Research, 115, 2010. 
O. Vasicek. An equilibrium characterization of the term structure. Journal of Financial Economics, 5(2): 

177 - 188, 1977. 

Q. H. Vuong. Likelihood ratio tests for model selection and non-nested hypotheses. Econometrica, 57: 
307-333, 1989. 



